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1.  Introduction 


The  maximin  algorithm  is  an  adaptive- array  algorithm  that  suppresses  interference  and 
thereby  supplements  the  inherent  processing  gain  of  a  spread-spectrum  system.  Originally 
developed  for  frequency-hopping  systems  (1),  the  algorithm  discriminates  between  the 
desired  signal  and  interference  on  the  basis  of  the  distinct  spectral  characteristics  of 
spread-spectrum  signals.  As  indicated  by  its  name,  the  maximin  algorithm  simultaneously 
maximizes  the  desired-signal  component  and  minimizes  the  interference  component  in  the 
despread  signal.  An  adaptive  array  using  the  maximin  algorithm  provides  a  direct-sequence 
system  with  a  high  degree  of  protection  against  strong  interference  that  would  overwhelm 
the  inherent  processing  gain.  The  main  advantages  of  the  maximin  algorithm  relative  to  its 
alternatives  (e.g.,  (2),  (3), (4))  are  that  it  does  not  require  training  sequences,  directional 
information,  decision- directed  adaptation,  or  elaborate  computations  such  as  eigenanalysis. 

The  basic  configuration  of  an  adaptive  array  for  spread-spectrum  systems  is  displayed  in 
figure  1.  The  spreading  waveform,  which  is  produced  by  a  synchronized  receiver,  enables 
the  despreading  of  each  copy  of  the  spread-spectrum  signal.  The  sample  values  of  the 
complex  envelopes  of  the  despread  signals  are  extracted  by  initial  processors  and  applied  to 
an  adaptive  processor  that  executes  the  maximin  algorithm.  In  a  direct-sequence  receiver, 
the  despread  desired  signal  has  a  narrowband  spectrum,  but  the  interference  has  a 
wideband  spectrum.  This  spectral  difference  is  exploited  by  the  maximin  algorithm  to 
estimate  the  interference  and  then  cancel  it. 

The  next  section  introduces  the  notation  and  provides  a  derivation  of  the  maximin  algorithm. 
In  section  3,  the  details  of  the  implementation  of  the  maximin  processor  in  a  direct-sequence 
system  are  explained.  The  convergence  analysis  of  section  4  establishes  bounds  on  the  adap¬ 
tation  constant  and  justifies  the  adaptation  sequence  used  in  the  maximin  algorithm.  Section 
5  presents  the  results  of  simulation  experiments.  The  derivation,  convergence  analysis,  and 
simulation  results  for  the  recursive  suppression  algorithm,  which  is  applicable  prior  to  code 
acquisition,  are  presented  in  section  6. 


2.  Derivation  of  Maximin  Algorithm 


The  desired  signal  and  the  interference  are  assumed  to  arrive  at  an  adaptive  array  of  N 
antennas.  The  desired  signal,  interference  signals,  and  thermal  noise  are  modeled  as 
independent  zero-mean,  wide-sense-stationary  stochastic  processes.  Let  x(i)  denote  the 
discrete-time  vector  of  the  complex  envelopes  of  the  N  antenna  outputs  after  each  one  has 
been  sampled,  despread,  and  filtered.  The  index  i  denotes  the  sample  number.  The  vector 
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To  demodulator 


Figure  1.  Architecture  of  adaptive  array  for  spread-spectrum  system.  SW  =  spreading 
waveform. 


x(i)  can  be  decomposed  as 


x(i)  =  s  (i)  +  n(z) 


(1) 


where  s(z)  is  the  vector  of  desired-signal  complex  envelopes,  and  n(/’)  is  the  vector  of 
interference  and  thermal-noise  complex  envelopes.  The  adaptive  filter  generates  a  weight 
vector  W  with  complex-valued  components.  The  output  of  the  adaptive  filter  is 


y(i)  =  WHx(i)  =  ya{i)  +  yn(i) 


where  H  denotes  the  complex  transpose  and 


3/,(i)  =  WHs(i), 


yn(i)  =  WHn(i). 


(2) 

(3) 


Let  E[  ]  denote  the  expected  value.  When  W  is  a  constant,  the  desired-signal  output 
power  is 

ps  =  £[|^(i)|2]=WHRsW  (4) 

where  Rs  is  the  desired- signal  correlation  matrix: 

Rs  =  E[s(i)sH{i)}.  (5) 


The  interference-plus-noise  output  power  is 


pn  =  WHRnW 


(6) 
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where  R„  is  the  interference-plus-noise  correlation  matrix: 


R  n  =  E[n(i)nH  (i)]. 

The  signal-to-interference-plus-noise  ratio  (SINR)  is 

_  ps  _  WhKsW 
P~Jn~  WHRnW  ' 

The  SINR  provides  the  performance  measure  that  the  adaptive  algorithm  seeks  to 
maximize. 


(7) 

(8) 


The  weight  vector  may  be  decomposed  as  W  =  Wr  +  jW /,  where  W /->  and  W r  are  the 
real  and  imaginary  parts  of  W,  respectively,  and  j  =  y/—l.  Let  V\vn  and  V w i  denote  the 
gradients  with  respect  to  Wr  and  W /,  respectively.  The  complex  gradient  is  defined  as 

=  -  (VifR  +  j'Vwi)  •  (9) 


The  maximin  algorithm  is  based  on  the  method  of  steepest  descent  (5),  (6).  In  this  method, 
the  weight  vector  is  changed  along  the  direction  of  the  negative  gradient  of  a  performance 
measure  that  is  to  be  minimized.  Separate  steepest-descent  equations  can  be  written  for 
W ft  and  W;.  Combining  these  equations  and  using  the  negative  SINR  as  the  performance 
measure,  we  obtain 

W(fc  +  1)  =  W(k)  +  p0(k)Vwp(k)  (10) 


where  k  denotes  the  weight  iteration  number,  /j0(/c)  is  a  scalar  sequence  that  controls  the 
rate  of  change  of  the  weight  vector,  and  V wp(k)  is  the  gradient  of  the  SINR  at  iteration  k. 
Equations  8  and  9  yield 


V  wp 


P 


RW 

Ps 


R»W 

Pn 


(11) 


Substitution  of  this  equation  into  equation  10  and  the  replacement  of  W  with  W (k)  gives 
the  steepest-descent  algorithm.  However,  this  algorithm  requires  knowledge  of  R  s  and  ps, 
which  are  difficult  to  estimate  directly.  To  obtain  a  more  practical  algorithm,  we  observe 
that  if  the  interference  and  noise  are  zero-mean  and  statistically  independent  of  the  desired 
signal,  then  the  input  correlation  matrix  is 

R;c  =  £’[x(i)xjH'(i)]  =  Rs  +  R„  (12) 


and  the  adaptive-filter  output  power  is 

Px  =  E[\y(i)  |2]  =  W^R,W  =ps  +  pn 


Substitution  of  equations  12  and  13  into  11  and  simplification  yields 


V  wp 


(p+i) 


R,W 

Px 


R„W 

Pn 


(13) 


(14) 
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which  may  be  approximated  by  observing  that  if  W  is  modeled  as  deterministic 
(nonrandom),  then  R^W  =  E[x.(i)y* (i)]  and  R„W  =  E[n(i)y*(i)\,  where  the  asterisk 
denotes  the  complex  conjugate.  Thus,  we  only  require  estimators  of  E[x.(i)y*(i)], 
E[n(i)yl(i)\,  px,  and  pn. 


A  further  simplification  that  ultimately  reduces  the  amount  of  computation  by  nearly  a 
factor  of  two  is  obtained  by  assuming  that  x(i)  is  obtained  by  filtering  samples  of  a 
continuous-time  vector  X(t)  of  complex  envelopes.  Each  component  of  X(t)  is  a  delayed 
version  of  the  first  component,  which  serves  as  a  reference  signal  and  may  be  expressed  as 

X1(t)  =  Xc(t)+jXs(t)  (15) 

where  Xc(t)  and  Xs(t)  are  the  real  and  imaginary  parts  of  Xi{t).  If  Xi(t)  is  modeled  as  a 
zero- mean,  wide-sense  stationary  process,  then  the  autocorrelations  of  Xc(t)  and  Xs(t)  are 
identical,  and  the  circular  symmetry  implies  that  E[Xc(t)Xs(t  +  r)]  =  —  E[Xs(t)Xc(t  +  r)], 
where  r  is  an  arbitrary  delay  (7),  (8).  From  these  relations,  it  follows  that 

-E[x(i)xT(i)]  =  0  (16) 

where  T  denotes  the  transpose.  The  adaptive-filter  output  can  be  decomposed  as 


y(i)  =  yr(i) +jyi(i) 


(17) 


where  yr(i)  and  yi{i)  are  the  real  and  imaginary  parts  of  y(i),  respectively.  If  the  weight 
vector  is  modeled  as  deterministic,  then  equations  2  and  16  imply  that 


xu 


ix"(s)W  +  lxT(i)W 


E[x(i)yr(i)\  =  E 

This  equation  and  equation  12  yield 

RiW  =  2  E[x.(i)yr(i)\. 


=  \e  [x(i)xff(i)]  W,  (18) 


(19) 


Similarly,  if  zero-mean,  wide-sense  stationary  interference  arriving  at  each  antenna  is  a 
delayed  version  of  the  interference  arriving  at  antenna  1,  and  if  the  zero-mean  thermal 
noise  is  independent  in  each  array  branch,  then 


£’[n(f)nT(f)]  =  0 


(20) 


and,  hence, 


R„W  =  2E[n(i)ynr(i)' 


(21) 


where  ynr(i)  is  the  real  part  of  yn{i).  Straightforward  calculations  using  equations  2,  16, 
and  13  yield  E[y2r(i)\  =  E[y‘f(i)\  and  E[yr(i)ys(i)\  =  0.  Therefore,  equations  13  and  17 
imply  that  px  =  2 E[yf(i)\.  Similarly,  pn  =  2 E[y^r(i)\.  Substitution  of  the  preceding  results 
into  sample  values  of  equation  14  yields 
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VM,p  —  (p  +  1) 


(22) 


E{x.yr }  _  E[ny  nr\  1 

Wr I  "  -Bfe]  J 

where  the  sample  index  has  been  omitted  for  simplicity. 


To  derive  the  maximin  algorithm,  let  px{k )  and  pn(k)  denote  estimates  of  E[y%\  and  E[y^r], 
respectively,  following  weight  iteration  k.  The  estimate  of  p  following  iteration  k  is  ’pik). 
Let  C x(k)  and  C n(k)  denote  estimates  following  iteration  k  of  the  input  correlation  vector 
E\xyr\  and  the  noise  correlation  vector  E[nynr],  respectively.  The  adaptation  sequence  is 
defined  as  a{k )  =  /i0(/c)[p(/c )  +  1].  Substituting  these  estimates  into  equations  22  and  10, 
we  obtain  the  basic  form  of  the  maximin  algorithm: 


W(fc  +  1)  =  W  (k)  +  a{k) 


C  x{k) 

_Px{k) 


C  njkY 

Pn{k )  _ 


k  >  0 


(23) 


where  W(0)  is  the  deterministic  initial  weight  vector.  As  the  adaptive  weights  converge, 
the  interference  components  of  Cx{k )  and  px(k )  decrease.  Thus,  the  first  term  within  the 
brackets  can  be  interpreted  as  a  signal  term  that  enables  the  algorithm  to  direct  the  array 
beam  toward  the  desired  signal.  The  second  term  within  the  brackets  is  a  noise  term  that 
enables  the  algorithm  to  null  interference  signals. 


The  adaptation  sequence  a(k )  should  be  chosen  so  that  I?[W(/c)]  converges  to  a  nearly 
optimal  steady-state  value.  It  is  also  intuitively  plausible  that  a{k )  should  decrease  rapidly 
as  if[W(/c)]  converges.  A  suitable  candidate  is 

Pn(k) 

t{k) 

where  t{k)  is  an  estimate  of  the  total  interference  and  noise  power  entering  the  array,  and 
a  is  the  adaptation  constant.  The  subsequent  convergence  analysis  and  simulation  results 
confirm  that  this  choice  is  effective  and  robust,  provided  that  the  adaptation  constant  is 
within  certain  numerical  bounds.  Another  adaptation  sequence  that  works  well  (1)  is 
a(k)  =  apn(k) / px(k) ,  but  the  best  choice  for  a  in  this  sequence  depends  on  the  SINR  at 
each  input  to  the  adaptive  filter. 

The  remaining  issue  is  the  choice  of  estimators  for  t(k),  C x(k),  C n(k),px(k),  and  pn{k).  The 
specific  nature  of  the  spread-spectrum  signals  allows  blind  estimates  to  be  made  without 
depending  on  known  steering  vectors  or  reference  signals. 


a(k)  =  a 


(24) 
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3.  Implementation  of  Processor 


The  maximin  processor  is  placed  after  the  despreading  subsystem  of  the  receiver.  For  the 
phase-shift  keying  (PSK)  modulation,  the  principal  components  of  the  initial  processor  in 
each  branch  following  an  antenna,  which  includes  the  despreading  subsystem,  are  shown  in 
figure  2.  The  front-end  devices  include  a  bandpass  filter  that  blocks  noise  outside  the  band 
occupied  by  the  spread-spectrum  signal.  Assuming  frequency  and  chip  synchronization  in 
the  receiver,  the  received  signal  is  downconverted  and  then  applied  to  a  chip  matched  filter. 
If  the  sampling  rate  equals  the  chip  rate,  the  desired  component  of  the  sampled  output  of 
the  chip  matched  filter  is  the  data-modulated  spreading  sequence.  If  code  synchronization 
of  the  spreading  sequence  has  been  established  in  the  receiver,  the  final  mixing  operation 
produces  the  branch  output  sequence  comprising  the  despread  desired  sequence,  spectrally 
spread  interference,  and  noise. 


Modulated 


Antenna 


Branch 

output 

sequence 


Figure  2.  Initial  processor  in  each  branch  of  the  receiver  for  a  direct-sequence  signal  with 
PSK. 


Although  frequency  synchronization  is  necessary  for  downconversion  of  the  received  signal 
to  baseband,  no  phase  synchronization  is  attempted  in  the  array  branches  because  relative 
phase  information  must  be  preserved  for  successful  beamforming  and  interference 
cancellation.  Phase  synchronization  is  ultimately  necessary  for  coherent  demodulation  of 
PSK.  The  random  phase  in  the  adaptive-filter  output  is  subsequently  removed  by  using  the 
output  of  a  carrier  synchronization  system. 

In  the  maximin  processor  shown  in  figure  3,  each  branch  output  sequence  is  applied  to  a 
baseband  filter  (BF), which  produces  a  component  of  x(i),  and  a  monitor  filter  (MF),  which 
estimates  the  interference.  For  an  adaptive  array  of  N  antennas,  the  outputs  of  N  pairs  of 
baseband  and  monitor  filters  are  applied  to  the  adaptive  filter,  the  output  of  which  is 
applied  to  a  digital  demodulator.  The  maximin  algorithm  seeks  to  maximize  the  SINR  at 
the  input  to  the  demodulator. 
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Branch  output 
sequences 


To 

demodulator 


Figure  3.  Maximin  processor. 

Figure  4  sketches  the  mainlobes  of  the  frequency  responses  of  the  baseband  and  monitor 
filters.  Each  baseband  filter  has  a  frequency  response  H(f)  with  one-sided  bandwidth 
B  «  1/TS,  where  Ts  is  the  data-symbol  duration.  The  N  baseband-filter  outputs  generate 
x(i).  Each  monitor  filter  has  a  frequency  response 

HAf)  =  \h{!  -  k)  +  \h(1  +  /„)  (25) 

where  the  center-frequency  offset  is  /o,  as  indicated  in  figure  4.  The  factor  1/2  ensures  an 
accurate  interference-power  estimate  when  the  interference  is  approximately  spectrally 
uniform  over  the  band  I/I  <fo  +  B. 


Magnitude 


Figure  4.  Frequency  responses  of  baseband  and  monitor  filters. 
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Each  baseband  filter  can  have  the  same  transfer  function  it  would  have  if  only  a  single 
antenna  were  used.  Thus,  since  the  despread  desired  signal  is  constant  during  a  symbol, 
this  signal  is  applied  to  a  matched  filter  with  a  rectangular  impulse  response.  The 
z-domain  transfer  function  of  the  matched  filter  is 

H(z)  =  l  +  z~1  +  ---  +  C  "  (26) 

-L  Z 

where  g  =  Ts/Tc  denotes  the  number  of  chips  per  symbol  or  the  processing  gain,  and  Tc  is 
the  chip  duration.  This  filter  has  a  null-to-null  bandwidth  equal  to  2 / gTc  =  2 jTs.  Each 
matched-filter  output  is  sampled  at  the  end  of  every  symbol  interval.  Let  Xi  (£)  denote  the 
vector  of  branch  output  sequences  applied  to  the  N  matched  filters.  The  vector  output  of 
these  filters  is 

i 

x(0  =  Xi(£)  (27) 

£=i—g+ 1 

where  l  is  the  index  of  the  chip-rate  input  samples  and  %  is  the  index  of  the  symbol-rate 
output  samples.  A  device  that  computes  this  sum  is  called  an  accumulator.  Similarly,  the 
monitor  filters,  each  of  which  has  frequency  response  are  bandpass  accumulators 

that  produce 

i 

n(i)  =  ^2  xi(£)cos(27 rf0Tc£)  (28) 

£=i—g+ 1 

which  is  the  vector  of  interference-plus-noise  estimates  used  to  generate  C n(k)  and  pn(k). 

The  despreading  of  the  direct-sequence  signal  spreads  the  spectral  density  of  the 
interference  over  the  entire  passband  of  the  monitor  filter  if  fo<(g  —  1  )/Ts.  Any  spillover 
or  spectral  splatter  of  the  desired-signal  spectrum  into  the  monitor  filter  may  lead  to  some 
degree  of  desired-signal  cancellation  by  the  adaptive  algorithm.  Thus,  /o  must  be 
sufficiently  large  to  prevent  significant  spectral  splatter. 


The  architecture  of  the  adaptive  filter  is  illustrated  in  figure  5.  The  vectors  applied  to  the 
adaptive  filter  are  x(i)  and  n(i),  and  m  symbol-rate  samples  are  taken  per  weight  iteration. 
The  adaptive  filter  produces  the  output 

yr(i)  =  -Re[WH(/c)x(i)],  i  —  km  +  1,  •  •  •  ,  (k  4*  l)m  (29) 

where  sample  %  is  taken  after  weight  iteration  k.  This  output  is  applied  to  the  demodulator 
and  is  used  in  the  estimators 

(k+l)m 

cx (k)  =  —  22  X(*)^(*)’  k  >  0  (30) 

i=km+ 1 


Figure  5.  Adaptive  filter  that  executes  maximin  algorithm. 


The  adaptive  filter  also  generates 

Vnrii)  =  i?e[WH(fc)n(i)],  i  =  km  +  1,  ■  •  •  ,  (k  +  l)m  .  (32) 

If  1)  /Ts  so  that  the  interference  power  is  spread  approximately  uniformly  over 

|/|  <  ,/o  +  B,  then  the  form  of  Hi (/)  indicates  that  the  interference  powers  at  the 
baseband-filter  and  monitor-filter  outputs  are  approximately  equal.  Since  the  noise 
components  are  zero-mean  and  independent  and  the  each  interference  component  is  a 
phase-shifted  version  of  the  component  in  a  reference  branch,  E[h(i)hH (i)]  &  £'[n(z)n'f/(i)]. 
Therefore, 

E[n(i)ynr(i)]  «  E[n(i)ynr(z)\  (33) 

and  a  suitable  estimator  of  the  interference-plus-noise  correlation  at  weight  iteration  k  is 

1  (k+l)m 

c n(k)  =  —  ^2  n(i)ynr(i),  k  >  0.  (34) 

i=km-\- 1 

Similarly,  a  suitable  estimator  proportional  to  the  interference-plus-noise  output  power  is 

1  (k+l)m 

Pn(k)  =  —  #nr(*)>  k  >  °-  (35) 

lit 

i=km-\- 1 


9 


An  estimator  of  the  total  interference-plus-noise  power  entering  the  array  is 

(k+l)m 

t(k)  =  —  ||  n(i)  ||2,  k  >  0.  (36) 

i=km+ 1 

where  ||  v  ||  denotes  the  Euclidean  norm  of  v.  These  estimators  complete  the  specification 
of  the  maximin  algorithm  given  by  equations  23  and  24. 

As  verified  by  simulation  experiments,  recursive  versions  of  the  preceding  estimators  merely 
slow  the  convergence  of  the  maximin  algorithm,  even  if  the  interference  statistics  are 
stationary.  Simulation  experiments  confirm  that  for  cyclostationary  spread-spectrum 
signals  and  tone  interference,  the  algorithm  simplification  stemming  from  the  assumption 
of  wide-sense  stationary  processes  causes  no  overall  performance  loss  other  than  a  slight 
slowing  of  convergence. 

In  a  multipath  environment,  separate  pairs  of  initial  and  maximin  processors  can  establish  a 
distinct  adaptive- array  pattern  for  each  resolved  path.  A  rake  combiner  can  then  maximize 
the  SINR  of  the  combined  signal  derived  from  all  these  patterns.  If  the  array  antennas 
are  sufficiently  close,  the  fading  across  the  array  is  completely  correlated,  and  the  adaptive 
processing  may  precede  the  channel  estimation  required  by  the  rake  combiner. 


4.  Convergence  Analysis 


Let  s{i)  denote  the  component  of  s(i)  derived  from  a  fixed  reference  antenna.  In  practical 
spread-spectrum  systems,  the  desired  signal  is  sufficiently  narrowband  that  its  copies  in  all 
the  branches  are  nearly  aligned  in  time.  Therefore, 


s  (i)  =  s(i)  S0  (37) 

where  So  is  a  steering  vector  of  complex  numbers  that  represent  the  relative  amplitudes 
and  phase  shifts  at  the  antenna  outputs.  The  substitution  of  this  equation  into  equation  5 
yields 

Rs  =  PsiSoSo  (38) 

where 

p„  =  B[|s(!)|'2].  (39) 

After  the  substitution  of  these  equations  into  equation  8,  the  maximization  of  the  SINR 
(appendix)  yields  the  optimal  weight  vector 

Wo  =  77R;1Sq  (40) 
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where  rj  is  an  arbitrary  constant.  The  maximum  SINR,  which  occurs  if  W  =  W0,  is 

p0=paiS*R?S0.  (41) 

The  highly  nonlinear  nature  of  the  maximin  algorithm  precludes  a  completely  rigorous 
convergence  analysis.  However,  with  enough  approximations,  the  convergence  of  the  mean 
weight  vector  to  W0  can  be  demonstrated  and  bounds  on  the  adaptation  constant  can  be 
derived.  We  assume  that  the  interference  is  wide-sense  stationary  and  m  is  large  enough 
that  equation  36  gives 

t(k)  ~  r  =  E[t(k )]  =  £'[||  n(i)  ||2]  ~  £'[||  n(i)  ||2]  =  tr(R,„)  (42) 

where  tr(  )  denotes  the  trace.  We  assume  that  after  a  number  of  algorithm  iterations  /c0, 
Px(k) /pn(k)  ~  p  +  1  ~  po  +  1.  Using  these  assumptions  in  equations  23  and  24,  the 
maximin  algorithm  is  approximated  by 

w (k  +  1)  =  W(fc)  +  -  -  cra(fc)1  ,  k>k0.  (43) 

r  Lido  +  ij 

Since  W (k)  does  not  depend  on  x(i)  and  n(i)  for  i  >  km  +  1,  we  make  the  approximation 
that  W (k)  is  statistically  independent  of  x(i)  and  n(i)  for  i  >  km  +  1.  We  obtain  from 
equations  30,  29,  16,  and  2  that 

E[Cx(k)]  =  E\x(i)yr(ij]  =  ^R^[W(fc)].  (44) 

Similarly,  equations  34,  33,  32,  20,  and  3  yield 

E[Cn(k)]  =  E[h(i)ynr(i)]  «  ^TLnE[W(k)}.  (45) 

Taking  the  expected  value  of  both  sides  of  equation  43,  substituting  equations  45,  44,  and 
12,  and  simplifying  algebraically,  we  obtain  the  approximate  recursive  equation  for  the 
mean  weight  vector: 

cy 

£[W(fc  +  l)]=  1-  D  B[W(t)],  fc>fc„  (46) 

where 

D  =  po Rn  -  Rs  =  po Rn  -  PsiS0So  .  (47) 

A  straightforward  calculation  using  equations  47  and  41  yields 

D  R^^o  =  0  (48) 
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which  indicates  that  R“  1 S0  is  an  eigenvector  of  D,  and  the  corresponding  eigenvalue  is  0. 
Since  D  is  Hermitian,  it  has  a  complete  set  of  TV  orthogonal  eigenvectors,  one  of  which  is 
R,n  1  So-  Since  p  <  po,  equation  8  implies  that  W/,RSW  <  poW^R^W.  Consequently, 
W/;D  W  >  0,  which  proves  that  D  is  positive  semidehnite  and,  hence,  has  TV  nonnegative 
eigenvalues.  Since  only  W  =  W0  gives  p  =  p0,  one  of  these  eigenvalues  is  zero,  and  the 
other  TV  —  1  eigenvalues  are  positive. 


To  solve  equation  46,  we  make  the  decomposition 

N 

E[W(k)]  =  p(fc)R-1S0  +  ^  (49) 

i=2 

where  each  a,i(k)  and  rj(k)  are  scalar  functions  and  e*  is  one  of  the  TV  —  1  eigenvectors 
orthogonal  to  R“1S0.  Substituting  this  equation  into  equation  46  and  using  the 
orthogonality  of  the  eigenvectors,  we  obtain 


&i(k  +  1)  — 


r){k  +  1)  =  r]{k)  =  rj(k0),  k  >  k0 

a\i 


1  - 


a,i(k),  2  <  i  <  TV,  k  >  kc 


(50) 

(51) 


2r(p0  +  1) 

where  A*  is  the  eigenvalue  corresponding  to  e*.  Assuming  that  rj(k0)  ^  0,  equations  40,  49, 
and  50  indicate  that  if[W(/c)]  — >  Wo  as  k  — >  oo  if  and  only  if  each  a,i(k )  — >  0.  The 
solution  to  equation  51  is 


di(k)  = 


1 


aXi 


k—ko 


di(ko),  2  <  i  <  TV,  k  >  k0. 


2r(p0  +  1) 

This  equation  indicates  that  eq(/c)  — >  0,  2  <  i  <  TV,  as  k  — »  oo  if  and  only  if 

ctA,- 


1  - 


2r(p0  +  1) 


<  1  ,  2  <  i  <  N. 


(52) 


(53) 


This  inequality  implies  that  the  necessary  and  sufficient  condition  for  the  convergence  of 
the  mean  weight  vector  is 

Mpo  +  1)  (54) 


0  <  a  < 


Xr 


where  Xmax  is  the  largest  eigenvalue  of  D. 


Since  the  sum  of  the  eigenvalues  of  a  square  matrix  is  equal  to  its  trace, 

N 

Xmax  A  E  A i  =  tr( D)  =  p0r  -  tr(Rs)  <  p0r.  (55) 

i= 1 

Substituting  this  bound  into  equation  54  and  simplifying  the  result,  we  obtain 

0  <  a  <  4  (56) 
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as  a  sufficient  (but  not  necessary)  condition  for  the  convergence  of  the  mean  weight  vector  to 
the  optimal  weight  vector.  Although  this  inequality  must  be  regarded  as  an  approximation 
because  of  the  approximations  used  in  its  derivation,  it  gives  at  least  rough  guidance  in  the 
selection  of  the  adaptation  constant.  The  fact  that  the  upper  bound  is  numerical  and  does 
not  depend  on  environmental  parameters  provides  support  for  the  choice  of  equation  24  as 
the  adaptation  sequence. 


5.  Simulation  of  Maximin  Algorithm 


In  the  simulation  experiments,  the  array  consists  of  four  omnidirectional  antennas  located 
at  the  vertices  of  a  square  or  in  a  uniform  linear  configuration.  Let  A  denote  the 
wavelength  corresponding  to  the  center  frequency  of  the  desired  signal,  which  is  3  GHz. 
The  edge  length  or  the  separation  between  adjacent  antennas  is  d  =  0.5A,  d  =  1.0A,  or 
d  =  1.5 A.  The  direct-sequence  signal  with  PSK  modulation  and  a  rectangular  chip 
waveform  arrives  from  a  direction  20  degrees  clockwise  from  the  perpendicular  to  one  of 
the  edges.  All  signals  are  assumed  to  arrive  as  plane  waves  with  no  fading.  The 
direct-sequence  signal  has  a  frequency  offset  equal  to  1  kHz  after  downconversion,  which 
models  imperfect  frequency  synchronization.  The  data-bit  and  spreading  sequences  are 
randomly  generated  for  each  simulation  trial  at  the  rates  of  100  kbps  and  10  Mbps, 
respectively.  The  processing  gain  is  20  dB.  Perfect  chip  and  spreading-sequence 
synchronization  in  the  receiver  are  assumed.  As  indicated  in  figure  2,  the  sampling  is  done 
once  per  spreading-sequence  chip.  The  thermal  noise  in  each  branch  output  is  modeled  as 
bandlimited  white  Gaussian  noise.  The  signal-to-noise  ratio  (SNR)  is  0  dB  at  the  input  of 
each  of  the  baseband  and  monitor  filters.  The  baseband  filters  are  accumulators.  The 
monitor  filters  are  bandpass  accumulators  offset  by  /o  =  400  kHz  to  prevent  contamination 
by  the  direct-sequence  signal.  The  maximin  algorithm  is  implemented  with  a  =  1.  A 
weight  iteration  occurs  after  each  ten  data  bits.  Each  simulation  trial  represents  the 
execution  of  the  algorithm  until  it  terminates.  For  each  trial,  the  initial  weight  vector  of 
the  adaptive  processor  is  W(0)  =  [1  0  0  0],  which  forms  an  omnidirectional  array 
pattern.  The  principal  system  parameter  values  are  listed  in  table  1. 

Each  of  1,  2,  or  3  interference  signals  is  a  tone  (continuous- wave  signal).  After  the 
downconversions,  the  tones  have  different  initial  phase  shifts  and  residual  frequency  offsets 
equal  to  10  kHz,  which  reflects  the  mismatch  of  the  tone  frequencies  and  the  carrier 
frequency  of  the  direct-sequence  signal.  Multiple  tones  do  not  add  coherently,  even  if  they 
have  the  same  carrier  frequencies,  because  they  arrive  from  different  directions  and  have 
different  initial  phase  shifts.  The  SINR  at  the  processor  output  is  calculated  after  each 
sample  time  and  then  averaged  over  all  samples  in  the  time  interval  between  a  weight 
iteration  and  a  preceding  weight  iteration  to  determine  the  SINR  at  each  weight  iteration. 


13 


Table  1.  Basic  system  parameters. 


Parameter 

Value 

Array  antennas 

4 

Array  geometry 

square  or  linear 

Center  frequency 

3  GHz 

Antenna  separation 

0.5A,  1.0A,  or  1.5A 

Direction  of  desired  signal 

20° 

SNR  in  each  branch  output 

0  dB 

Modulation 

PSK 

Data  rate 

100  kbs 

Filter  band  widths 

100  kHz 

Monitor  filter  offset 

400  kHz 

Frequency  offset  of  desired  signal 

1  kHz 

Chip  and  sampling  rate 

10  Mbs 

Bits  per  weight  iteration 

10 

Weight  iterations  per  trial 

100 

Adaptation  constant 

a  =  1 

Initial  weight  vector 

[1  0  0  0] 

Frequency  offset  of  interference 

10  kHz 

The  SINR  is  observed  to  fluctuate,  but  tends  to  gradually  increase  until  it  reaches  a 
steady-state  condition  with  a  smaller  residual  fluctuation. 

Figure  6  illustrates  the  SINR  variation  versus  the  weight  iteration  number  for  a  typical 
simulation  trial  in  which  one  interference  tone  arrives  at  a  50°  angle  with  an 
interference-to-signal  ratio  (ISR)  equal  to  30  dB.  Let  0  denote  an  arrival  angle  defined  as 
the  angle  in  the  counterclockwise  direction  from  the  normal  to  one  of  the  array  edges.  Let 
S  (9)  denote  the  array  response  vector ,  which  is  the  array  response  to  an  ideal  plane  wave 
arriving  at  angle  9  (5),  (8).  The  array  gain  pattern  after  weight  iteration  k  is 

G(6,k)  =\WH(k)S(9)\2 .  (57) 

The  array  gain  pattern  at  the  end  of  the  simulation  trial  of  figure  6  is  depicted  in  figure  7. 
A  null  deeper  than  —20  dB  in  the  direction  of  the  interference  signal  and  a  mainlobe 
slightly  displaced  from  the  direction  of  the  desired  signal  have  formed  along  with  other 
grating  nulls  and  grating  lobes. 

The  results  of  15  representative  simulation  experiments  using  the  parameter  values  of  table 
1  are  summarized  in  table  2.  Each  experiment  comprises  50  trials  with  100  weight 
iterations  per  trial.  The  first  column  gives  the  arrival  angles  of  1,  2,  or  3  interference 
signals  relative  to  the  desired-signal  direction.  The  ISR  for  each  interference  signal  is 
10  dB.  The  second  column  gives  the  array  type:  square  with  d  =  1.0A,  square  with 
d  =  1.5A,  or  linear  with  d  =  0.5A.  The  SINRs  expressed  in  decibels  for  the  last  20  weight 
iterations  of  all  the  trials  are  used  to  compute  the  steady-state  SINR  and  the  standard 
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Figure  6.  SINR  variation  in  typical  simulation  trial  for  maximin  algorithm,  one  interfer¬ 
ence  tone  at  50°,  and  ISR  =  30  dB. 


Figure  7.  Array  gain  pattern  at  end  of  typical  simulation  trial  for  maximin  algorithm, 
one  interference  tone  at  50°,  and  ISR  =  30  dB. 
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Table  2.  Simulation  results  for  interference  tones,  each  with  ISR  =  10  dB. 


Arrival  angles 
of  interf.  (°) 

Array 

Type 

Steady-state 
SINR  (dB) 

Standard 
Dev.  (dB) 

Crossing 

Number 

60 

Sq,  1.0A 

25.81 

1.49 

3.2 

60,  -40 

Sq,  1.0A 

22.84 

2.22 

8.9 

60,  85 

Sq,  1.0A 

24.53 

1.60 

5.0 

60,  -40,  85 

Sq,  1.0A 

24.39 

1.71 

5.5 

30 

Sq,  1.0A 

20.10 

1.57 

16.1 

60 

Sq,  1.5A 

25.13 

1.46 

4.2 

60,  -40 

Sq,  1.5A 

25.38 

1.50 

3.9 

60,  85 

Sq,  1.5A 

24.40 

1.47 

5.4 

60,  -40,  85 

Sq,  1.5A 

24.26 

1.52 

5.4 

30 

Sq,  1.5A 

22.95 

1.99 

7.1 

60 

Linear 

25.75 

1.43 

3.3 

60,  -40 

Linear 

24.98 

1.46 

4.1 

60,  85 

Linear 

24.81 

1.52 

4.4 

60,  -40,  85 

Linear 

24.67 

1.59 

4.8 

30 

Linear 

20.25 

1.55 

10.1 

deviation  of  the  SINR.  The  final  column  gives  the  crossing  number ,  which  is  the  average 
number  of  weight  iterations  required  for  the  SINR  to  exceed  20  dB.  The  crossing  number 
provides  a  rough  measure  of  the  relative  time  required  for  convergence  to  the  steady  state. 

Table  2  indicates  that  beamforming  and  interference  cancellation  are  achieved  in  a  wide 
variety  of  scenarios.  Additional  interference  signals  do  not  necessarily  slow  the  algorithm 
convergence  or  lower  the  steady-state  SINR.  The  reason  is  that  as  the  array  forms  one 
pattern  lobe  or  null,  it  tends  to  form  additional  grating  lobes  and  nulls  and  thereby  may 
counteract  the  maximin  algorithm.  The  results  for  an  interference  signal  with  a  30°  arrival 
angle  and  a  square  array  illustrate  the  limitations  imposed  by  the  resolution  of  the  array. 
The  resolution ,  which  is  the  angular  separation  between  the  interference  and  desired-signal 
sources  that  can  be  accommodated  without  a  large  performance  degradation,  decreases  as 
the  array  aperture  increases.  The  square  array  with  d  =  1.0 A  and  the  uniform  linear  array 
with  d  =  0.5A,  despite  having  an  array  aperture  equal  to  1.5A,  have  insufficient  resolution 
for  this  interference.  Increasing  the  antenna  separation  in  a  square  array  to  1.5 A  provides  a 
large  improvement  in  the  convergence  speed  and  the  steady-state  SINR. 

Table  3  lists  the  simulation  results  for  a  square  array  with  d  =  1.0A,  one  interference  tone 
with  a  60°  arrival  angle,  and  a  variable  ISR.  Each  experiment  comprises  50  trials  with  100 
weight  iterations  per  trial.  The  parameters  of  table  1  are  used  except  that  there  are  15  bits 
per  weight  iteration  when  ISR  =  35  dB.  The  enormous  interference-cancellation  capability 
of  the  maximin  processor,  which  supplements  the  inherent  processing  gain  of  the 
direct-sequence  signal,  is  evident.  As  the  ISR  increases  from  10  dB  to  35  dB,  the 
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Table  3.  Simulation  results  for  a  square  array  with  d  =  1.0 A  and  one  interference  tone 
with  a  60°  arrival  angle. 


ISR  (dB) 

Steady-state 
SINR  (dB) 

Standard 
dev.  (dB) 

Crossing 

number 

10 

25.81 

1.49 

3.2 

20 

24.22 

1.44 

6.2 

30 

23.64 

1.44 

18.2 

35 

23.66 

2.61 

44.3 

convergence  slows.  The  delay  corresponding  to  the  crossing  number  is  1.82  ms  when  ISR  = 
30  dB,  but  partially  because  of  the  decrease  in  the  weight  iteration  rate,  this  delay 
increases  to  6.645  ms  when  ISR  =  35  dB. 

As  the  ISR  increases  beyond  30  dB,  it  is  found  that  the  convergence  of  the  SINR  to  the 
steady  state  for  all  trials  requires  an  increase  in  the  number  of  bits  per  weight  iteration.  In 
general,  an  increase  in  the  number  of  bits  per  weight  iteration  enhances  stable  convergence 
and  reduces  the  standard  deviation  of  the  SINR.  However,  the  steady-state  SINR  is  reduced, 
and  the  convergence  is  slower.  The  choice  of  10  bits  per  weight  iteration  in  table  1  represents 
a  reasonable  compromise  in  the  performance  objectives  for  ISRs  that  do  not  exceed  30  dB. 


6.  Acquisition 


A  direct-sequence  signal  cannot  be  despread  until  the  receiver  can  establish  a  local 
spreading  sequence  synchronized  with  the  received  spreading  sequence  (8).  Prior  to  initial 
sequence  synchronization  or  code  acquisition,  the  usual  adaptive-array  algorithm  for 
interference  suppression,  such  as  the  maximin  algorithm,  might  be  ineffective  against 
strong  interference  because  the  relative  strengthening  of  the  desired  signal  by  means  of 
despreading  is  unavailable.  The  wide  variety  of  techniques  that  require  only  a  single 
antenna  output  to  reject  narrowband  interference  (8)  rely  on  sequence  synchronization 
and,  thus,  are  not  useful  prior  to  code  acquisition.  Adaptive  beamforming  potentially 
would  increase  the  desired-signal  power  relative  to  the  interference,  but  the 
direction-of-arrrival  estimation  and  the  avoidance  or  rejection  of  beams  in  the  directions  of 
strong  interference  signals  might  be  excessively  time  consuming. 

Typical  receivers  for  direct-sequence  signals  only  activate  a  single  antenna  during  code 
acquisition  (2).  A  proposed  algorithm  to  enhance  acquisition  by  using  the  full  array  (9)  is 
extremely  computationally  intensive  and  requires  short  spreading  sequences,  which  are 
usually  unacceptable.  A  more  recent  proposed  algorithm  (10)  that  uses  the  normalized 
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least-mean-square  algorithm  does  not  have  these  disadvantages,  but  neither  this  algorithm 
nor  the  other  one  can  prevent  beamforming  in  the  direction  of  strong  interference.  Thus, 
the  algorithms  may  amplify  the  interference  instead  of  cancelling  it.  Another  approach  is 
to  use  the  recursive  suppression  algorithm  prior  to  code  acquisition  for  the  purpose  of 
interference  suppression  without  significant  desired-signal  cancellation.  The  reduced 
interference  in  the  adaptive-processor  output,  which  is  applied  to  the  acquisition  system, 
facilitates  the  code  acquisition  process. 

The  recursive  suppression  algorithm  is  designed  to  cancel  all  received  signals.  If  the 
interference  signals  are  stronger  than  the  desired  signal,  they  are  cancelled  first,  which 
enables  the  acquisition  system  to  achieve  code  synchronization.  The  algorithm  terminates 
before  it  has  enough  time  to  cancel  the  desired  signal.  If  acquisition  has  been  verified,  the 
receiver  activates  both  the  tracking  system  and  the  maximin  algorithm.  If  acquisition  has 
not  been  verified,  the  recursive  suppression  algorithm  is  repeated  periodically  with  an 
initial  omnidirectional  weight  vector.  If  the  desired  signal  is  stronger  than  the  interference, 
then  the  algorithm  is  not  useful,  but  the  acquisition  system  only  needs  to  process  the 
output  of  a  single  antenna,  which  it  does  in  the  interims  between  algorithm  repetitions. 

Each  antenna  output  is  applied  to  a  separate  initial  processor  that  produces  a  branch 
output  sequence,  as  illustrated  in  figure  2  for  a  direct-sequence  signal  with  PSK 
modulation.  The  first  mixer  downconverts  the  received  bandpass  signal  to  a 
complex-valued  baseband  signal.  After  the  downconversion,  a  chip  matched  filter  and 
sampler  produce  the  data-modulated  spreading  sequence.  This  sequence  is  mixed  with  the 
local  spreading  sequence  to  produce  a  chip-rate  branch  output  sequence  that  becomes  the 
despread  sequence  after  code  acquisition  has  occurred. 

As  illustrated  in  figure  8  for  complex-valued  signals,  branch  output  sequences  are  applied 
directly  to  an  adaptive  filter  that  executes  the  recursive  suppression  algorithm  and  is 
embedded  within  a  serial-search  or  digital  matched-filter  acquisition  correlator  (8).  The 
adaptive  filter  applies  the  weights  computed  by  the  recursive  suppression  algorithm  to  the 
branch  output  sequences.  The  complex- valued  adaptive-filter  output  is  applied  to  an 
accumulator  in  a  serial-search  system  or  digital  transversal  filters  in  a  matched-filter 
system.  The  acquisition  correlator  produces  a  decision  statistic  that  is  compared  to  a 
threshold.  Because  of  the  mixing  with  the  local  spreading  sequence,  both  the  noise  and  the 
interference  at  the  input  to  the  accumulator  or  transversal  filters  approximate  white 
Gaussian  processes,  and  the  analytical  results  of  (8)  for  calculating  the  acquisition 
performance  statistics  are  applicable. 
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Figure  8.  Adaptive  filter  and  acquisition  correlator. 


6.1  Derivation  of  Recursive  Suppression  Algorithm 

The  recursive  suppression  algorithm  can  be  derived  by  the  method  of  steepest  descent  with 
the  constrained  minimization  of  the  output  power  of  an  adaptive  filter  as  the  performance 
criterion.  Let  x(z)  denote  the  discrete-time  input  vector,  each  component  of  which  is  a 
branch  output  sequence.  If  x(i)  is  modeled  as  a  wide-sense  stationary  process  and  the 
weight  vector  W  produced  by  the  adaptive  filter  is  modeled  as  deterministic,  then  the 
adaptive-filter  output  is 


y(i)  =  WHx(i) 

(58) 

and  the  mean  output  power  is 

Px  =  W^R,  W 

(59) 

where  the  input  correlation  matrix  is 

Ra:  =  12[x(i)xH(i)]. 

(60) 

The  constraint  that  the  norm  of  the  weight  vector  must  remain  constant  ensures  that  the 
optimal  weight  vector  is  not  the  zero  vector.  For  mathematical  tractability,  the 
performance  criterion  is  prescribed  as  the  minimization  of  the  mean  output  power  Px 
subject  to  the  quadratic  constraint 


W  ||2=  WHW  =  1. 


(61) 


To  derive  the  optimal  weight  vector  by  the  method  of  Lagrange  multipliers,  we  minimize 
the  real  scalar 

H  =  P,-1(\\  W||2-l)  (62) 

where  7  is  a  scalar  Lagrange  multiplier.  Substituting  equation  59,  we  find  that  the  complex 
gradient  of  H  is 

WwH  =  2RXW  -  27W.  (63) 
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A  necessary  condition  for  the  minimum  is  determined  by  setting  this  equation  equal  to  the 
zero  vector,  which  yields  the  following  equation  that  must  be  satisfied  by  the  optimal 
weight  vector  W0: 

R,Wo  =  7W0.  (64) 

Thus,  W0  must  be  a  unit-norm  eigenvector  of  R,  with  eigenvalue  7.  Substituting 
W  =  Wo,  equations  64,  and  61  into  59,  we  obtain  Px  =  7.  Since  R„  is  Hermitian  and 
nonnegative  definite,  7  >  0  and  Px  >  0.  Therefore,  Px  is  minimized  when  7  =  Xrrun ,  the 
smallest  eigenvalue  R , .  We  conclude  that  the  optimal  weight  vector  is  a  unit-norm 
eigenvector  of  R.c  with  eigenvalue  A min- 


The  method  of  steepest  descent  for  discrete-time  systems  and  equation  63  give  the 
recursive  equation  for  the  unconstrained  weight  vector: 

W'(fc  +  1)  =  W'(fc)  -y0(k)WwH 

=  [1  +  27/i0(/c)]W/(/c)  -  2y0(k)KxW'(k)  (65) 


where  Ho(k)  is  a  real  sequence  that  regulates  the  convergence.  Dehning 

fii(k)  =  2/i0(/c)/[l  +  2 7/i0(/c)],  the  recursive  equation  for  the  constrained  weight  vector 

W(lfe)  =  W \k)/  ||  W '(k)  ||  is 


w ,,  |  W(fc)-/x1(fc)RxW(fc) 
1  J  ||  W{k)  -  fj,1{k)RxW(k) 


(66) 


To  convert  this  equation  into  a  practical  algorithm,  we  observe  that  if  W(/c)  is  modeled  as 
deterministic,  then  RxW(/c)  =  £,[x(f)xiJ(f)]W(/c)  =  E[x.(i)y*(i)],  where  i  is  the  sample 
index,  k  is  the  weight  iteration  index,  and  the  asterisk  denotes  the  complex  conjugate. 
Thus,  we  only  require  an  estimator  for  E,[x(i)|/*(i)]. 


A  further  simplification  that  ultimately  reduces  the  amount  of  computation  by  nearly  a 
factor  of  two  is  obtained  by  assuming  that  x(i)  is  obtained  by  sampling  a  continuous-time 
vector  X(t)  of  complex  envelopes.  Each  component  of  X(t)  is  a  delayed  version  of  the  first 
component,  which  is  modeled  as  a  zero- mean,  wide-sense  stationary  process.  As  shown  in 
section  2,  this  model  implies  equation  16  and  that 


RX.W  =  2 E[x(i)yr(i)  . 


(67) 


Substituting  this  equation  into  equation  66  and  dehning  a{k)  =  2pi(/c),  we  obtain  the 
recursive  suppression  algorithm-. 


W(fc  +  1) 


W(k)-a(k)Cx(k) 

W(k)-a(k)Cx(k) 


(68) 


where  C x(k)  is  an  estimate  of  E[x(i)yr(i)]  and  W(0)  is  the  specified  unit-norm  initial 
weight  vector. 
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The  vector  applied  to  the  adaptive  filter  is  x(i),  and  m  chip-rate  samples  are  taken  per 
weight  iteration.  The  adaptive  filter  produces  the  output 

(69) 


yr(i)  =  Re[WH  (k)x.(i)],  i  =  km  +  !,•••  ,  (k  +  l)m 


where  sample  i  is  taken  after  weight  iteration  k.  This  output  is  used  in  the  estimator 

1  (fc+l)m 

c x(k)  =  —  ^2  x(%r(*)>  k>  0.  (70) 

771 

i=km+l 


The  adaptation  sequence  a(k)  is  chosen  so  that  a{k)Cx(k)  has  dimensionless  components 
and  so  that  the  mean  weight  vector  converges.  As  shown  subsequently,  a  suitable  choice  is 


a(k) 


a 

tx{k ) 


(71) 


where  a  is  the  adaptation  constant,  and  an  estimator  of  the  total  power  entering  the 
adaptive  filter  is 

(k+l)m 

tx(k )  =  —  ^2  II  x(®)  II2’  k  -  0  ’  (72) 

i=km+ 1 

Simulation  experiments  confirm  that  for  cyclostationary  spread-spectrum  signals  and  tone 
interference,  the  recursive  suppression  algorithm  suffers  no  degradation  from  the  simplifica¬ 
tion  stemming  from  the  assumption  that  X\[t)  is  wide-sense  stationary.  As  has  been  verified 
by  simulation,  recursive  versions  of  equations  70  and  72  merely  slow  the  convergence  of  the 
algorithm,  even  if  the  signal  statistics  are  stationary. 


6.2  Convergence  of  Mean  Weight  Vector 


To  analyze  the  convergence  of  the  mean  weight  vector  of  the  recursive  suppression 
algorithm,  the  estimate  tx(k)  is  approximated  by  its  mean  value  E[tx(k)\  =  fr(R.,.).  The 
norm  of  V(/c)  =  W [k)  —  a{k)Cx(k )  in  equation  68  is  approximated  by  ||  E[V(k)]  ||.  Both 
of  these  approximations  become  more  accurate  as  k  increases.  Since  W (k)  does  not  depend 
on  x(i)  for  i  >  km  +  1,  we  make  the  approximation  that  W (k)  is  statistically  independent 
of  x(i)  for  i  >  km  +  1.  Therefore,  equation  16  implies  that 

E[Cx(k)]  =  E[x.{i)yr(i)}  «  ^E[x(i)x.H(i)}E[W(k)]  =  ^KxE[W(k)}.  (73) 


Taking  the  expected  values  of  both  sides  of  equation  68  and  using  the  approximations,  we 
obtain 


a;w(a-  +  i); 


(I  -  CR,)£[W(fc)] 
(I  -  (Rx)E[W(k)] 


k>  0 


(74) 
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where 


a 


(75) 


2//'!  Jt  ,:)  ' 

By  mathematical  induction,  equation  74  implies  that 


E[W(k)] 


(i  -  CR,)fcw(o) 
(i  -  CR,)fcw(o) 


k  >  o. 


(76) 


Since  R  ,  is  Hermitian,  it  has  a  complete  set  of  orthonormal  eigenvectors.  For  an  TV  x  N 
matrix  R,. ,  let  Sm]n  denote  the  space  spanned  by  the  N  —  M  orthonormal  eigenvectors 
associated  with  A min,  the  smallest  eigenvalue,  and  let  Sm  denote  the  space  spanned  by  the 
M  remaining  orthonormal  eigenvectors.  An  arbitrary  initial  weight  vector  can  be 
decomposed  as 

M 

W(0)  =  u  +  ^Av,  (77) 

i=  1 

where  u  is  in  Smin ,  the  {/3i}  are  scalars,  and  v*,  i  =  1,  2,  •  •  •  ,  M,  is  the  orthonormal  set  of 
eigenvectors  in  Sm-  The  probability  that  an  arbitrarily  chosen  nonzero  W(0)  will  have 
u  =  0  is  zero,  so  we  ignore  this  possibility.  Since 


UHV<  =  vf Vj  =  0,  i,  j  =  1,2,  -  -  -  ,M, 
the  substitution  of  equations  77  into  76  yields 

(1  -  CAmm)fcU  +  Yh=1  A( 1  -  C Ai)^* 


£?[W(fc  +  1)  =  - 

(l-CA™„)2fe  ||uP+EfiA2(l-CA*)2fc 

where  A,:  is  the  eigenvalue  associated  with  v*. 


-,1/2 


(78) 


(79) 


If  every  coefficient  of  v*  is  eventually  exceeded  in  magnitude  by  the  coefficient  of  u  as  k 
increases,  then  E[W(k  +  1)]  converges  to  an  optimal  value.  Thus,  if 

|l-CAi|  <  |l-CAmin|,  i  =  1,2,  ■  ■  ■  ,  M  (80) 

then 

lirn  E[W(k)]  =  .  (81) 

k—>  oo  U 

Conversely,  if  equation  81  holds  for  arbitrary  W(0)  and,  hence,  an  arbitrary  choice  of  the 
Pi,  i  =  1,  2,  •  •  •  ,  M ,  then  equation  80  must  be  satisfied.  Thus,  equation  80  is  a  necessary 
and  sufficient  condition  for  the  convergence  of  the  mean  weight  vector  to  an  optimal  value. 


Since  R„  is  nonnegative  definite,  A m%n  and  each  A*  are  nonnegative.  By  separately 
evaluating  equation  80  for  the  cases  in  which  <  1  and  (X *  >  1,  we  verify  that  a 
sufficient  convergence  condition  is 

o  <  C  <  t - (82) 

^ max  '  A min 
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where  Xmax  is  the  largest  eigenvalue  of  Rx .  Since  Xmax  +  Xmin  <  tr( R;c),  it  follows  from 
equations  82  and  75  that 

0  <  a  <  4  (83) 

is  a  sufficient  condition  for  the  convergence  of  the  mean  weight  vector  to  an  optimal  value. 
The  fact  that  the  bounds  are  numerical  and  do  not  depend  on  environmental  parameters 
provides  support  for  the  choice  of  equation  71  as  the  adaptation  sequence.  Simulation 
experiments  indicate  that  the  convergence  rate  tends  to  increase  as  a  increases  within  the 
bounds  of  equation  83. 

6.3  Simulation  of  Recursive  Suppression  Algorithm 

Simulation  experiments  were  conducted  to  discover  the  extent  to  which  interference  can  be 
cancelled  while  avoiding  the  cancellation  of  the  desired  signal.  In  most  of  the  simulation 
experiments,  the  array  consists  of  four  omnidirectional  antennas  located  at  the  vertices  of  a 
square.  The  edge  length  is  equal  to  the  wavelength  A  corresponding  to  the  center  frequency 
of  the  desired  signal,  which  is  3  GHz.  The  direct-sequence  signal  with  PSK  modulation  and 
a  rectangular  chip  waveform  arrives  from  a  direction  20  degrees  clockwise  from  the 
perpendicular  to  one  of  the  edges.  All  signals  are  assumed  to  arrive  as  plane  waves  with  no 
fading.  The  direct-sequence  signal  has  a  frequency  offset  equal  to  1  kHz  after 
downconversion,  which  models  imperfect  frequency  synchronization.  The  signal  samples 
are  multiplied  by  the  local  spreading  sequence,  which  is  not  synchronized  with  the 
spreading  sequence  of  the  received  direct-sequence  signal.  However,  the  minor  effects  of  the 
imperfect  sampling  times  and  the  self- interference  (8)  are  not  modeled  in  the  simulation  for 
simplicity.  The  multiplication  and  the  presence  of  data  modulation  produce  a 
data-sequence  component  of  the  branch  output  sequence  that  is  independent  of  the 
spreading  sequence  imposed  on  the  interference  and  noise  by  the  multiplication.  Both 
sequences  are  randomly  generated  for  each  simulation  trial  at  the  chip  rate  of  10  Mbps.  As 
indicated  in  figure  2,  the  sampling  is  done  once  per  spreading-sequence  chip.  The  thermal 
noise  in  each  branch  output  is  modeled  as  bandlimited  white  Gaussian  noise.  The  SNR  is 
assumed  to  be  only  0  dB  because  of  the  absence  of  successful  despreading  of  the  desired 
signal  and  data-bandwidth  filtering  prior  to  code  acquisition,  which  also  makes  the 
processing  gain  irrelevant.  Weight  iterations  occur  once  per  100  spreading-sequence  chips. 
Each  simulation  trial  represents  the  execution  of  the  algorithm  until  it  terminates.  For 
each  trial,  the  initial  weight  vector  of  the  adaptive  processor  is  W(0)  =  [1  0  0  0],  which 
forms  an  omnidirectional  array  pattern.  The  principal  system  parameter  values  are  listed 
in  table  4. 

The  adaptive-array  output  is  assumed  to  be  applied  to  a  serial-search  acquisition  system 
designed  to  achieve  code  acquisition  within  400  /jls  with  a  very  high  probability  if  its  input 
SINR  is  maintained  above  —6  dB  (cf.  (8),  chapter  4).  Thus,  the  number  of  weight 
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Table  4.  Basic  system  parameters. 


Parameter 

Value 

Array  antennas 

4  at  vertices  of  square 

Center  frequency 

3  GHz 

Array  edge  length 

wavelength 

Direction  of  desired  signal 

20° 

SNR  in  each  branch  output 

0  dB 

Modulation 

PSK 

Frequency  offset  of  desired  signal 

1  kHz 

Sampling  and  chip  rate 

10  Msamples/s 

Chips  per  weight  iteration 

100 

Weight  iterations  per  trial 

100 

Adaptation  constant 

a  =  0.08 

Initial  weight  vector 

[1  0  0  0] 

Frequency  offset  of  interference 

10  kHz 

iterations  per  trial  is  100,  which  corresponds  to  an  algorithm  duration  T0  =  1  ms.  This 
duration  allows  600  [is  for  enough  interference  suppression  to  attain  an  SINR  above  —6  dB 
and  another  400  /is  for  acquisition  to  be  completed.  The  choice  of  the  adaptation  constant 
as  a  =  0.08  is  found  to  produce  an  algorithm  convergence  rate  that  provides  rapid 
interference  cancellation  with  little  desired-signal  cancellation  when  the  initial  SINR  is 
sufficiently  low.  Simulation  experiments  indicate  that  the  adaptation  constant  should  be 
chosen  to  be  inversely  proportional  to  the  algorithm  duration.  Thus,  if  Tq  were  equal  to  4 
ms,  then  a  =  0.02  would  be  a  good  choice. 

Each  of  1,  2,  or  3  interference  signals  is  a  tone  (continuous- wave  signal),  which  generally 
degrades  a  direct-sequence  system  more  than  noise  interference  with  the  same  power  (8). 
After  the  downconversions,  the  tones  have  different  initial  phase  shifts  and  residual 
frequency  offsets  equal  to  10  kHz,  which  reflects  the  mismatch  of  the  tone  frequencies  and 
the  carrier  frequency  of  the  direct-sequence  signal.  The  SINR  at  the  processor  output  is 
calculated  after  each  sample  time  and  then  averaged  over  all  samples  in  the  time  interval 
between  a  weight  iteration  and  a  preceding  weight  iteration  to  determine  the  SINR  at  each 
weight  iteration. 

Figure  9  illustrates  the  SINR  variation  versus  the  weight  iteration  number  for  a  typical 
simulation  trial  in  which  one  interference  tone  arrives  at  a  60°  angle  with  an  ISR  of  10  dB. 
For  this  trial,  the  SINR  is  above  —6  dB  for  more  than  800  fis,  making  acquisition  almost 
assured.  The  array  gain  pattern  after  weight  iteration  k  is  defined  by  equation  57.  Figure 
10  depicts  the  gain  pattern  at  the  end  of  the  simulation  trial  of  figure  9.  A  null  deeper  than 
—20  dB  in  the  direction  of  the  interference  source  has  formed  along  with  a  grating  null  and 
a  partial  grating  null. 
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Figure  9.  SINR  variation  in  typical  simulation  trial  for  recursive  suppression  algorithm, 
one  interference  tone  at  60°,  and  ISR  =  10  dB. 


Figure  10.  SINR  variation  in  typical  simulation  trial  for  recursive  suppression  algorithm, 
one  interference  tone  at  60  °,  and  ISR  =  10  dB. 
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Table  5.  Simulation  results  for  interference  tones,  each  with  ISR  =  10  dB,  and  square 
array  with  d  =  X. 


Arrival  angles  of 
Interference  (°) 

Final 

SINR  (dB) 

Standard 
dev.  (dB) 

Relative 
gain  (dB) 

60 

-0.68 

0.47 

9.73 

60,  -40 

-3.22 

1.70 

16.82 

60,  85 

-3.26 

1.23 

16.78 

60,  -40,  85 

-3.76 

2.22 

26.24 

60,  -40,  105 

-1.46 

1.22 

29.54 

40 

-5.39 

0.56 

5.02 

The  results  of  6  representative  simulation  experiments  using  the  parameter  values  of  table 
4  are  summarized  in  table  5.  Each  experiment  comprises  50  trials  with  100  weight 
iterations  per  trial.  The  first  column  gives  the  arrival  angles  of  1,  2,  or  3  interference 
signals  relative  to  the  desired-signal  direction.  The  ISR  for  each  interference  signal  is 
10  dB.  The  SINRs  expressed  in  decibels  for  the  last  40  weight  iterations  of  all  the  trials  are 
used  to  compute  the  final  SINR  and  the  standard  deviation  of  the  SINR,  which  are  listed 
in  the  second  and  third  columns.  The  fourth  column  lists  the  gain  of  the  final  SINR 
compared  with  the  initial  SINR. 

The  first  row  of  the  table  shows  the  system  response  to  a  single  interference  signal 
separated  by  40°  from  the  desired  signal.  As  illustrated  in  figures  9  and  10,  the  interference 
signal  is  largely  suppressed  while  the  desired  signal  endures  much  less  cancellation.  The 
second  and  third  rows  indicate  that  there  is  some  decrease  in  the  final  SINR  due  to  the 
presence  of  a  second  equal-power  interference  signal,  but  both  interference  signals  are 
substantially  cancelled.  The  fourth  and  fifth  rows  illustrate  the  surprising  result  that  the 
final  SINR  and  the  standard  deviation  (SD)  do  not  necessarily  degrade  when  a  third 
equal-power  interference  signal  is  present  because  both  grating  lobes  and  grating  nulls  are 
naturally  generated  as  interference  signals  are  nulled.  The  large  values  of  SD  in  the  second, 
third,  and  fourth  rows  indicate  that  for  some  trials  the  SINR  is  not  maintained  above 
—6  dB  for  400  /is,  as  desired.  Therefore,  algorithm  repetitions  may  be  needed  for  code 
acquisition  to  be  achieved. 

The  final  row  illustrates  the  limitations  imposed  by  the  resolution  of  the  array,  which  is  the 
minimum  angular  separation  between  the  interference  and  desired-signal  sources  that  can 
be  accommodated  without  significant  performance  degradation.  If  the  array  edge  length  or 
antenna  separation  is  increased  to  d  =  1.5 A,  the  resolution  improves,  and  the  final  SINR 
increases  by  3.78  dB  relative  to  the  SINR  for  d  =  A.  However,  for  the  scenario  of  the  first 
row,  the  increased  antenna  separation  causes  the  formation  of  a  partial  grating  null  near 
the  direct-sequence  signal  that  reduces  the  final  SINR  by  3.51  dB. 
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If  interference  arrives  at  angles  in  the  interval  (—90°,  +90°),  then  the  resolution  can  be 
improved  without  deep  partial  gating  nulls  in  the  interval  by  using  a  uniform  linear  array 
of  4  antennas  with  d  =  0.5A.  The  improved  resolution  is  due  to  the  aperture  of  the  linear 
array,  which  is  1.5A.  Table  6  lists  the  results  for  simulation  examples  similar  to  those  in 
table  5  except  for  the  array  geometry. 


Table  6.  Simulation  results  for  interference  tones,  each  with  ISR  =  10  dB,  and  uniform 
linear  array  with  d  =  0.5A. 


Arrival  angles  of 
Interference  (°) 

Final 

SINR  (dB) 

Standard 
dev.  (dB) 

Relative 
gain  (dB) 

60 

-1.17 

0.47 

9.24 

05 

o 

1 

o 

-0.44 

1.64 

19.60 

60,  85 

-1.05 

1.38 

19.00 

60,  -40,  85 

-0.77 

1.35 

29.23 

40 

-0.01 

0.48 

10.40 

Table  7  lists  the  simulation  results  for  the  parameter  values  of  table  4,  one  interference 
tone  with  a  60°  arrival  angle,  and  a  variable  ISR.  Each  experiment  comprises  50  trials  with 
100  weight  iterations  per  trial.  As  the  ISR  increases  from  3  dB  to  30  dB,  the  convergence 
of  the  recursive  suppression  algorithm  slows,  but  the  final  SINR  and  the  relative  gain 
increase  monotonically  while  the  standard  deviation  remains  small  until  ISR>  20  dB. 
Although  an  SINR  loss  or  an  acquisition  failure  may  occur  when  ISR<  6  dB,  the  receiver 
can  acquire  the  direct-sequence  signal  by  processing  the  output  of  a  single  antenna  between 
algorithm  repetitions  or  whenever  a  low  ISR  is  known  to  exist. 


Table  7.  Simulation  results  for  one  interference  tone  at  60°  and  square  array  with  d  —  A. 


ISR  (dB) 

Final 

SINR  (dB) 

Standard 
dev.  (dB) 

Relative 
gain  (dB) 

3 

-4.97 

1.05 

-0.20 

6 

-2.67 

0.69 

4.32 

10 

-0.68 

0.47 

9.73 

20 

-0.02 

0.99 

20.02 

30 

-3.84 

2.85 

26.16 

The  results  in  table  7  are  for  a  fixed  adaptation  constant  a  =  0.08,  which  reflects  the  fact 
that  the  interference  environment  is  generally  unknown  and  time  varying.  However,  the  best 
choice  for  the  adaptation  constant  depends  on  the  ISR,  and  a  =  0.08  is  close  to  optimal  only 
for  ISR  ~  10  dB.  The  low  final  SINR  and  high  SD  for  ISR  =  30  dB  are  due  to  the  insufficiently 
rapid  cancellation  of  this  strong  an  interference  signal.  If  the  adaptation  constant  is  chosen 
to  be  a  =  0.16  when  ISR=  30  dB,  then  the  final  SINR  =  1.31  dB,  SD  =  0.44  dB,  and  the 
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relative  gain  is  31.31  dB.  At  the  other  extreme,  if  ISR  =  3  dB  and  a  =  0.04,  then  the  final 
SINR  =  —3.74  dB,  SD  =  0.47  dB,  and  the  relative  gain  is  1.03  dB.  These  results  indicate 
the  benefit  of  slowing  the  algorithm  convergence  when  the  interference  is  not  much  stronger 
than  the  desired  signal. 


7.  Conclusions 


An  adaptive  array  using  the  maximin  algorithm  provides  a  direct-sequence  system  with  a 
high  degree  of  protection  against  strong  interference  that  would  overwhelm  the  inherent 
processing  gain.  The  main  advantages  of  the  maximin  algorithm,  relative  to  its 
alternatives,  are  that  it  does  not  require  training  sequences,  decision-directed  adaptation, 
directional  information,  or  elaborate  computations. 

The  substantial  interference-suppression  capability  of  a  direct-sequence  system  depends  on 
its  ability  to  first  achieve  code  acquisition.  However,  in  the  presence  of  strong  directional 
interference,  perhaps  from  several  sources,  the  code-acquisition  process  may  be  severely 
impaired  or  completely  thwarted.  The  recursive  suppression  algorithm  appears  to  be  the 
only  existing  adaptive-array  algorithm  that  can  often  suppress  strong  interference  in  a  direct- 
sequence  system  enough  to  enable  rapid  code  acquisition.  The  algorithm  can  be  executed 
periodically  over  short  time  intervals  to  provide  repeated  acquisition  opportunities.  During 
the  interims  between  algorithm  executions,  the  output  of  a  single  antenna  can  be  processed 
to  enable  acquisition  when  the  interference  is  weak  or  absent. 
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A.  Appendix  -  Optimization  of  the  SINR 


The  signal-to-interference-plus-noise  ratio  (SINR)  at  the  adaptive-filter  output  is 

_W"RSW 

Po  W"R,„W  '  ^ 

Since  their  definitions  ensure  that  Rs  and  R„  are  Hermitian  and  nonnegative  dehnite, 
these  N  x  N  matrices  have  complete  sets  of  orthonormal  eigenvectors  with  nonnegative 
real-valued  eigenvalues.  The  noise  power  is  assumed  to  be  positive,  which  implies  that  R„ 
is  positive  dehnite  and  has  positive  eigenvalues.  The  spectral  theorem  of  linear  algebra 
indicates  that  R„  can  be  expressed  as 


N 

Rn  ^  ^ 

i= 1 


(A-2) 


where  A*  is  an  eigenvalue  and  e,  is  the  associated  eigenvector. 


To  derive  the  weight  vector  that  maximizes  the  SINR  with  no  restriction  on  R  s ,  we  dehne 
the  Hermitian  matrix 

L 


A  ^  '/\ieiei  . 

i=  1 

(A-3) 

Direct  calculations  verify  that 

R„  =  A2 

(A-4) 

and  the  inverse  of  A  is 

L 

(A-5) 

The  matrix  of  A  specifies  an  invertible  transformation  of  W  into  the  vector 

V  =  AW. 

(A-6) 

We  dehne  the  Hermitian  matrix 

C  =  A-1RSA-1. 

(A-7) 

Then  equations  A-l,  A-4,  A-6,  and  A-7  indicate  that  the  SINR  can  be  expressed  as 

a 

Rayleigh  quotient : 

VHC  V 

(A-8) 

P°  ~  Y  2  • 

Let  £i,  ■  ■  ■  An  and  u1;  u2,  •  •  • 

,  uN  denote  the  eigenvalues  and  corresponding 

orthonormal  eigenvectors  of  C. 

If  V  is  expanded  as 

N 

V  =  ^;Ui  (A-9) 

i— 1 
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where  the  {fej}  are  coefficients,  then 


N  N 

VHC  V  =  ^  \bi\%  <  imax  Y,  M2  =  tmaX  ||  V  II2  (A-10) 

i=  1  i=  1 

where  imax  is  the  largest  eigenvalue.  Therefore,  po  <  Imax ■  Direct  substitution  indicates 
that  po  is  maximized  by  V  =  7711,  where  u  is  an  eigenvector  of  C  associated  with  eigenvalue 
Imax,  and  77  is  an  arbitrary  constant.  Thus,  the  maximum  value  of  po  is 


Pmax  ^ max  • 


(A-ll) 


From  equation  A-6  with  V 
the  SINR  is 


7711,  it  follows  that  an  optimal  weight  vector  that  maximizes 


W0  =  //A  xu  . 


(A-12) 


The  purpose  of  an  adaptive-array  algorithm  is  to  adjust  the  weight  vector  to  converge  to 
the  optimal  value,  which  is  given  by  equation  A-12  when  the  maximization  of  the  SINR  is 
the  performance  criterion. 


If  a  desired  signal  with  input  power  psi  is  sufficiently  narrowband,  then  equation  38 
indicates  that  Hs=psiSoSQ  ,  where  So  is  the  steering  vector.  Substitution  into  equation  A- 7 
yields 

C  =  psl ,FFh  (A-13) 

where 

F  =  A_1S0  .  (A- 14) 

The  factorization  explicitly  shows  that  C  is  a  rank-one  matrix,  which  has  only  one  nonzero 
eigenvalue.  By  direct  substitution,  it  is  found  that  the  eigenvector  associated  with  the 
nonzero  eigenvalue  is 

u  =  F  =  A_1S0  (A-15) 

and  the  eigenvalue  is 

tmax=Psi  ||  F  ||2  .  (A-16) 

Substituting  equations  A-15  into  A-12  and  then  observing  that  R”1  =  A~2,  we  obtain  the 
optimal  weight  vector: 

W0  =  vR-'So  (A-17) 

where  77  is  an  arbitrary  constant.  The  maximum  value  of  the  SINR,  obtained  from 
equations  A-ll,  A-16,  and  A-14  is 


Pr, 


=  p.s,;S"  R„  1  Sc 


(A-18) 
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